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' We review properties of gapless states. We construct a model where a stable 

breached pair (gapless) state is realized. 

o ■ 
o 

' 1. Introduction 

Qh! Motivated by recent experiments in cold atoms ^ and by questions in QCD 
^ • at high densities ^, we consider here superfluid fermion systems, in partic- 
le . ular exotic superfluids. As suggested by Liu and Wilczek, exotic phases of 
'"^ ' matter in superfluid fermion systems involve coexistence of normal Fermi 

liquid and superfluid components Superfluid properties are described by 

^ I a nonzero condensate, A = {T{^^x)^^x'))) ^ 0, being the order parame- 

' ter, while a single quasiparticle dispersion crosses the momentum axis (free 



Fermi surface) leading to a gapless mode; thus this phase is sometimes called 
a "gapless" superconductive phase. Due to nonzero Fermi condensate the 
ground state is a superfluid in a classical sense {i.e. it has zero viscosity). 
On a microscopic scale, one can envision a momentum separation in exotic 
superfluids. For species with noticeably different Fermi momenta. Cooper 
pairing takes place around the Fermi surfaces, but there is no pairing in 
the momentum region between surfaces (the breach). Thus in this work we 
use the term "breached pair" (BP) superfluidity. 

To obtain exotic superfluids we consider pairing between two different 
fermion species whose Fermi surfaces do not match. This possibility arises 
in several situations: (1) Spin-up spin-down electrons in an ordinary su- 
perconductor placed in a uniform magnetic field undergo Zeeman splitting, 
leading to a mismatch in Fermi momentum. As found by Sarma in the 
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1960's ^, an exotic superconducting ground state should arise at large mo- 
mentum mismatch when iibH > A. Before that, however, the first-order 
phase transition from the superconducting state to the normal state takes 
place at hbH = A/\/2. Placing a superconductor in a spacially varying 
magnetic field or adding paramagnetic impurities with a strong spin-fiip 
electron-impurity scattering amplitude stabilizes the gapless superconduc- 
tor ^. (2) Recent experiments in cold atomic fermion gases trapped in an 
optical lattice and operating near Feshbach resonance deal with a mixture 
of two hyperfine spin components of alkali atoms. By changing the scat- 
tering length one can go from the regime of Bose- Einstein condensation to 
BCS superfluidity, which is of interest in this work. Laser lattice involves 
counterpropagating laser beams, that together generate a standing light 
wave leading to different AC Stark shifts for the spin-up and spin-down 
components. Using methods of 'engineering' various lattice systems and by 
tuning effective masses one can produce exotic phases ®. (3) In strongly 
interacting quark matter at high baryon densities and low temperatures, 
different flavors of quarks pair and form color superconductors ^. Here a 
mismatch in Fermi momenta arises due to a nonzero strange quark mass, 
nis ^ 0, and is triggered by imposing a charge neutrality condition. At 
intermediate density (2 3 nuclear densities), an exotic state may arise 
which links the CFL and nuclear matter phases. 

Plan: First, we give a general analysis of gapless states. Then, we show 
how to realize a stable BP superfluid state. 

2. Gapless state and its stability 

We consider the mean-flield analysis of a model with two species of fermions 
A, and B of differing masses niA < tub and in the absence of inter- 
action with different Fermi momenta pA < Pb- The Hamiltonian is 
H = J cPp/ {2tt)^ V'ApV'Ap + Sp'^^pi'Bj^ + Hi with attractive interac- 
tion Hi = -g J d^p/{27Tfd^q/{2n)^^\^ijl,_pi;Bgi^A-q {g > 0) and = 
p'^/2mA -HA,£p^ p^/2mB - Mb, where pa = y/2rnjjl^, pb = V2tobMb- 
At the mean-field level, the condensate A = / d^p/ {2Tr)^ {^pB pi' A- p) is a 
c-number, which permits to diagonalize the Hamiltonian. As a result, the 

quasiparticle excitations are = s~ ± Sp + with = (e^ ±5^ )/2; 
they contain mixture of A-particle and B-hole excitations. We minimize the 
thermodynamic potential Q = H — ^AnA — IJ'SnB, i.e. dQ/dA = 0, to find 
the gap parameter A. There are two non-trivial solutions. First one with 
larger gap corresponds to a fully gapped BCS state, where have oppo- 
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site signes for all momenta. Though ^ ns at 5 = 0, in the presence of 
interaction particles redestribute so that the occupation numbers become 

equal ua = fiB = n with n = 1/2 ^1 — Sp /'\Js:t'^ + A^^ ; the BCS con- 
densation energy is the largest for pA = Pb ■ Second solution, obtained by 
Sarma in 60's, has smaller gap. Interaction is not strong enough to pull both 
Fermi surfaces together, leaving a breach region with single occupancy by B- 
particles. Pairing and superfluidity takes place primarily around the smaller 
Fermi surface, while there is normal component for momenta in the breach 
region where have the same signs (separation in momentum). Points 
where = 0, i.e., pf ^ = +p|)/2 ± l/2^{pl - p\)^ - 16mAmBA\ 
give the free Fermi surfaces leading to gapless modes and to free fermion 
liquid component. Occupation numbers are fiA = fiB = n for < p < pi 
and p2 < p, and n_B = 1 ua = for pi < p < P2, Figure 1. 

We consider nontrivial superfluid solutions of the gap equation in grand 
canonical and canonical ensembles. 

I. Grand canonical ensemble. 

We fix the chemical potentials /i^, Mb and minimize the thermodynamic 
potential Q,s = H — HAnA — ^J-BnB over all ground state superconducting 
wave functions, mm(^5|0|^'5). We obtain, apart from the trivial solution 
A = 0, two nontrivial solutions: the fully gapped BCS and the gapless 
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Figure 2. Solutions of tlie gap equation and corresponding pressures as functions of the 
Fermi momenta mismatch at fixed chemical potentials (left), and gap equation solutions 
and energies at fixed particle numbers (right). 



Sarma solution ^. As a function of the Fermi momenta mismatch 5p = 
{Pb - Pi)/4VmIrrri, they are Assc = Aq and l^sarma/ = \/'^x - 1 
with X = Sp/Ao, Figure 2. Since at a given 6p and g, Abcs > Asarma, 
BCS state wins, i.e., system prefers the BCS ground state over the Sarma 
state. At fixed Fermi momenta pa,Pb, the thermodynamic potential as a 
function of the gap has two minima - normal (N) and the BCS (absolute 
min) states, and one maximum - the Sarma state. Hence the Sarma state 
is mctastablc. 

We consider the pressure versus the Fermi momentum mismatch for 
states which are solutions of the gap equation, where normalized pres- 
sure is defined through the condensation free energy as Ps = — ((fis) — 
(r2o))/|((^^BCs) — (^^o))| where (fio) is the free energy of the normal state at 
5p = Q and the BCS condensation energy is {^bcs) — i^o) = — -^(0)Ao/2, 
Ao is the BCS gap and N{0) is the density of states at the Fermi sur- 
face. In the leading order A ^ Sp <^ pa,Pb (i-C, when all quantities arc 
written as expansions near the Fermi surface), pressures are Pbcs = 1 for 
< a; < 1, Psarma = '^x^ " (1 " 2x)2 for 1/2 < X < 1 and Pat = 2x'^ for 
a; > with x = Sp/Ao, Figure 2. State with maximum pressure (or min- 
imum condensate energy) wins. For pA = Pb there is the BCS state. As 
we add B-particles, we create stress, and pressure of the BCS state Pbcs 
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drops relatively to the pressure of the normal unpaired state Pn ~ Sp'^. 
When Pbcs ^ Pn < 0, there is a first order phase transition at A 7^ 
from the superconducting to the normal state. The BCS superconductor 
is destroyed when win from condensation energy besomes less than loss in 
energy needed to pull the Fermi surfaces together to create the BCS. Sarma 
state is tangent to normal state at A = 0, and has always lower pressure 
than normal state; hence Sarma state is unstable. 
II. Canonical ensemble. 

We fix the particle numbers allowing the chemical poten- 

tials to change, and minimize the Helmholtz energy over all possible 
superconducting ground states with constraint of fixed rij i = A,B, 
min{'i' slHl'i' s)ni-const ^- At a single point when ua — ns there is the 
BCS state with I^bcs = Aq. At tia 7^ there is Sarma state with 
decreasing gap bs 5n = ub — tla increases, Asarma/ = a/1 — 2a; where 
X = Sp/ Ao ^ 6n. The energy of Sarma state is lower than the energy 
of the normal state, Esarma < -E-at; thus Sarma state is stable. Imposing 
neutrality condition a stable gapless superconducting state was obtained in 
the QCD context by Shovkovy et. al. ^ . 

We came to different conclusions about stability of Sarma state using 
grand canonical and canonical ensembles. This difference in stability analy- 
sis can be resolved by considering mixed phase, which is a mixture in space 
of two (or more) homogeneous states. Bedaque et. al. suggested to con- 
sider a mixture of the BCS and normal states, separated in x-space. They 
found that the energy of the mixed state is lower than the energy of Sarma 
state, Emixed < Esarma', thus Sarma state is unstable with respect to decay 
into a mixed state. We confirm their findings. We define the normalized 
condensation energy as Es = {{Hs) — {Hn)) /\{{Hbcs) — {Hn))\ where 
{Hbcs) - {Hn} = -iV(0)A2/2. In the leading order A ~ (5p < Pa,Pb, 
condensation energies are Ebcs = —1 for a; = 0, Esarma = ~(1 ~ 2a;)^ for 
< a; < 1/2, Ernixed = —(1 — V^x)'^ and E^ = 0, x = Sp/Ao and Aq is 
the BCS gap, Figure 2. 

Allowing mixed states in our ansatz of trial ground state wave functions 
^s, metastable Sarma state decays (rolls down) into a mixture of the BCS 
and normal states. Figure 3. (In a mixed state, pressures and chemical 
potentials of composite states are equal, hence there are two equal minima 
^Bcs = ^n)- Generally, it is difficult to include mixed states in variational 
ansatz in the grand canonical ensemble. 

We conclude, that in grand canonical and canonical ensembles Sarma 
branch of the gapless solutions is unstable. 
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Figure 3. Thermodynamic potential as a function of the gap, and positions of each 
state. Sarma state decays into a mixed BCS and Normal states depicted on the right 
panel. 

In we show that conclusion about stabihty of a state is the same 
in any ensemble used. In particular, there is one-to-one correspondence 
between a state at fixed particle number(s) and the state that minimizes 
the thermodynamic potential O in grand canonical ensemble. Thus, there 
is always a stable state in the grand canonical ensemble that satsfics the 
constraint. Imposing constraint (over particle numbers) cannot stabilize 
the system. Practical guide is to look for a stable (gapless) solution in the 
grand canonical ensemble. 

3. Breached paired superfluid state for a finite-range 
interaction 

Our goal is to construct a stable breached paired state which has coexisting 
superfluid and gapless components. We use an idea that existance of the 
gapless modes depend on the momentum structure of the gap A(p). There 
should be two distinct regions in momentum space: first one where Ap is 
large enough to support the superfluid, and second one where Ap is small 
enough that pairing does not appreciably affect the normal free-fermion 
behavior. To garantee stability, the phase must also have higher pressure 
than the normal state. We realize such BP states in two examples 
I. Cut-off interaction. 

We impose a cut-off interaction such that it supports the BCS-like 

pairing for p < pA and it allows free dispersion relations for p > pa, 
accomodating the excess of B-particles and leading to gapless modes. 
We construct this state by minimizing the thermodynamic potential, 
min{^s\H — fiA^A — A^b'^bI^'s), where H includes the cut-off interac- 



tion -g J £p/(27rfd^q/{27rff{p)f{q)i:\jjl i;B-gtPAg with f{p) = 1 



for p > PA and f{p) = for p < pa- The gap parameter, deflned 
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Figure 4. Gap, dispersions and occupation numbers for the cut-oflF (left) and the two- 
body potential (right) interactions. 

as A = gj (fp/{2Tr)^f{p){^Bp'4'A-p), satisfies the gap equation A = 

l/2g J d^q/ (27r)^ A/(g) / y^ep ^ + A^ where momentum integration is per- 
formed outside the breach region. The occupation numbers nA,nB show 
the evidance that it is a breached paired state, Figure 4. This state is an 
absolute minimum of the thermodynamic potential, hence we obtained a 
stable BP state. 

II. Spherically symmetric static two-body potential. 

With attractive potential V{x — x'), interaction is Hj = 
J (fip/{2-K)^(fiq/{2'KYV{p — q}4^ApVB-pVB-qil'Aq and the gap parameter 
acquires a momentum dependence, Ap — J d?q/{2'KYV{p ~ q){'>pBq^A-q) ■ 

The gap equation is written Ap = 1/2 / (fq/{2T^Yv{p - q)Aq/^ej ^ + A^, 

and quasiparticle dispersion relations are E^=e-± ^Je+'^ + Al- Wc take 
a gaussian potential for numerical simulations. Due to the BCS instability. 
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Ap picks at the effective Fermi surface given by the pole of the gap equation 
at A = 0, = 0. Therefore Ap supports the BCS-hke pairing around po, 
and allows free dispersion relations, and hence free Fermi surfaces, outside 
the breached region, Figure 4. It is, however, difficult to varify that this 
state is an absohite minimum of the thermodynamic potential since instead 
of a number. A, we have a function Ap in the variational ansatz. 

We performed minimization of the thermodynamic potential numeri- 
cally using different potentials. Generally, there is a central strip of fully 
gapped BCS phase about pA = Pb with normal unpaired phase outside. De- 
pending on the parameters of interaction, these phases may be separated 
by a region of gapless BP superfluid phase. Figure 5. 

Conditions to have BP phase are as follows. At pA ~ Pb there is 
standard BCS, which is a stable fully gapped solution. By adjusting the 
chemical potentials so as to increase the Fermi surface ps, we stress the 
system and low the pressurerelative to the normal phase. Eventually, either 
before or after a transition to a BP state, the pressure becomes negative 
and there is a first order phase transition to the normal phase. At the 
point just before transition: if Ap^ is sufficiently large, the state is fully 
gapped (BCS) and no BP state will occur; if Ap^^ is small, then it will 
not appreciably affect the dispersions and one finds a gapless Fermi surface 
coexisting with the superfluid phase. As long as Ap falls off sufficiently 
quickly, one can choose large ratio ms/mA 1 so that the transition will 
occur with Ap^ small enough to support the BP phase. States shown at 
Figure 5 have uiB/'mA = 10. For a wider in q-space interaction, larger 
mass ratio is needed. 




Figure 5. Phase diagram of possible homogeneous phases in coordinates of the Fermi 
momenta (pa^Pb) for the cut-off (left) and two-body potential (right) interactions. 
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4. Conclusion 

We considered a Fermi system with weak attractive interaction between 
species A and B, where the BCS state forms at equal number densities, 
riA = nB- What is the ground state of this system when ua «b? There 
is a range of parameters, where a breached pair phase exist. Breached pair 
state is a homogeneous phase where supcrfluid and normal components 
coexist. This state is stable and can be found provided there is a momentum 
structure of interaction and large enough mass ratio of two species. 
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